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1 Introduction 



The quantum two-time Green function technique has provided a very powerful 
tool in condensed matter physics for exploring the equilibrium and transport 
properties of a wide variety of many-body systems. Within this framework, the 
equation of motion method (EMM) and the spectral density method (SDM) 
allow to obtain reliable approximations to treat typically unperturbative prob- 
lems [1,2,3]. The pioneering introduction of the two-time Green functions and 
the EMM in classical statistical mechanics by Bogoljubov and Sadovnikov [4], 
has opened the concrete possibility to describe classical and quantum systems 
on the same footing. Nextly, a classical version of the SDM (CSDM) has been 
also formulated in extensive [5,6,7] and nonextensive [8] classical statistical 
mechanics and applied to classical magnetic chains [5,6,7,8]. The last method, 
which seems to present many advantages with respect to the most conven- 
tional EMM [3,5,6,7,8], offers a robust instrument for systematic and well 
tested approximations [3] to explore the macroscopic properties of classical 
many-body systems as well. Unfortunately, both the mentioned methods have 
not received the due consideration in the classical context and further develop- 
ments and applications are desiderable. Along this direction, in this paper we 
apply the CSDM to investigate the thermodynamic properties of a highly non- 
trivial d- dimensional classical spin-S* Heisenberg ferromagnet with long-range 
interactions decaying as r~ p {p > d) with the distance r between spins in the 
presence of an external magnetic field. The choice of this model, besides its in- 
trinsic theoretical interest, is motivated also by its relevance in the description 
of many materials [9,10,11] of experimental and technological importance and, 
in particular, in the understanding of the critical behavior of magnetic systems 
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when long-range exchange interactions are believed to play an important role. 
In general, realistic microscopic magnetic models with long-range interactions 
are difficult to be studied in a reliable way. Since the exact solution for the 
spin-1/2 quantum Heisenberg chain with an inverse- square exchange (p = 2) 
and without an external field was found independently by Haldane [12] and 
Shastry [13], only a limited lot of information has been acquired about the 
low-temperature properties, the existence of long-range order (LRO) and the 
related critical behavior for different values of the decaying exponent p and 
the dimensionality d. 

Recently, the critical properties at finite temperature of <i-dimensional quan- 
tum Heisenberg models, with interactions of the type here considered (here 
often named "long-range spin models") have been studied using microscopic 
techniques [14,15,16,17] and Monte Carlo simulations [18]. A summary of the 
known features sounds as follows. The one- and two-dimensional long-range 
quantum spin-1/2 Heisenberg ferromagnets in absence of an external magnetic 
field were investigated by Nakano and Takahashi using the so called modified 
spin- wave theory [14] and the Schwinger-boson mean-field approximation [15]. 
Further information were derived for the ci-dimensional case by means of the 
EMM for the two-time Green functions using the Tyablikov decoupling proce- 
dure [16]. Monte Carlo simulations for the two-dimensional quantum spin-1/2 
Heisenberg model have been also performed for 2 < p < 6 [18]. This scenario, 
has been recently enriched by an extension [17] of the Mermin- Wagner theo- 
rem [19] for the existence of ferromagnetic (FM) LRO at finite temperature 
in quantum Heisenberg and XY models in d(— 1, 2) dimensions with r~ p - and 
oscillatory- interactions. 

Classical long-range spin-s Heisenberg FM models have attracted great atten- 
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tion, too. It has been proved that LRO exist in d(— 1,2) dimension when 
d < p < 2d [20,21] and is destroyed at all finite temperatures for p > 2d 
[22,23,24]. Similar results were obtained for the spherical model [25] and the 
present scenario of the critical properties is largely based on renormalization 
group calculations for the classical n- vector model [26,27]. Classical long-range 
antiferromagnetic (AFM) models have been studied less extensively. The avail- 
able rigorous results [20,21,22] suggest orientational disorder at all finite tem- 
peratures when p > 2d, but no theorem exists entailing existence or absence 
of LRO for d < p < 2d. Monte Carlo simulations have been also performed 
for both d(= 1, 2) -dimensional classical FM (for p = 2d [28]) and AFM (for 
p = 3/2 and p = 3 with d = 1 and d = 2 respectively [29]) Heisenberg 
long-range models. The results confirm that FM-LRO survives at finite tem- 
perature provided d < p < 2d and allow to conjecture that no AFM-LRO 
exists at all finite temperatures for p > d. Spin-wave studies [29] agree with 
last conjecture but no definitive statement can be drawn at the present stage. 
In a quite recent work [30] the thermodynamics and the critical properties of 
the classical long-range spin-S* Heisenberg FM chain in the presence of an ex- 
ternal field have been systematically studied using the CSDM to lowest order 
of approximation. The results are in good agreement with previous available 
analytical and numerical investigations. As mentioned before, here we use the 
same method to explore a d-dimensional Heisenberg FM model. As we shall 
see, a rich phase diagram appears where all the above mentioned scenario is 
reproduced and extended. 

The paper is organized as follows. In Sec. 2 we introduce the model and the 
moment equations (ME's) for the spectral density (SD) to lowest order of 
approximation. The main low-temperature properties are studied in Sec. 3 
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and the existence of LRO in different regions of the (d — p)-plane is also shown 
in a transparent way. Sec. 4 is devoted to the critical properties and the low- 
temperature paramagnetic susceptibility of the model. Finally, in Sec. 5, some 
concluding remarks are drawn. 

2 The model and the moment equations for the spectral density 

A classical rf-dimensional spin-5 1 Heisenberg ferromagnet with long-range in- 
teractions is described by the Hamiltonian: 

i N N 

H = -~Y. Ji&- Sj-hJ2 S ?- (!) 

z ij=i i=i 

Here, N is the number of sites of the hypercubic lattice with unitary spacing, 
{Si]i = 1, ...,N} are the classical spins, h is the external magnetic field and 
the spin-spin interaction, in view of the thermodynamic limit as iV — > oo, 
is assumed to be = J/rfj, where = |r» — r^l and J > measures 
the strength of the coupling. The extreme case p — > oo corresponds to the 
standard nearest-neighbor interaction while the mean field approximation is 
obtained when p = (replacing J by J/N). For this type of interaction, the 
thermodynamical limit iV — > oo is well defined only for p > d, while for p < d 
the ground state of the system has an infinite energy per particle as iV — > oo 
and the conventional statistical mechanics cannot be directly applied. 

The classical spins model (1) can be appropriately described by the set of 2N 
canonical variables {ipj, 5?} where (fj is the angle between the projection of the 
spin vector Sj in the (x — y)-plane and the x axis. For practical calculations 
we find convenient to introduce the new variables Sf — ± iSj, so that 
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Sj = (Sj) 2 + SfSj . Then, the Hamiltonian (1) can be rewritten as 

H = E J ( k ) (^ 5 -k + SJP'k) " (2) 

k 



involving the Fourier components of the spins and the exchange interaction 
defined by: 



N N 



S k = E ' ' k , 'Sr J (k) = E e*- (ri " Pi) ^, (3) 

where k denotes a wave vector in the ci-dimensional Fourier space. The sum 
in Eq. (2) is restricted to the first Brillouen zone (1BZ) of the lattice. The 
Poisson brackets for the spin Fourier components relevant for us, are: 

{'S'kS'Sk'} = i^k+k'; {^k'^k'} = ~^^k+k' ■ (4) 



In the context of the CSDM [5,6,7,8,30], all the thermodynamical properties 
of the model can be derived using the SD defined by: 

A k (u;) = -i ({SZ k , (t)}) u = -i dte™ ({Sl k , S£ (*)}) , (5) 

where the classical dynamical variables depend on time through the conju- 
gate canonical coordinates {</?,• (t) , S? (t) — ...ivj. Here A (t) = e lLt A (0), 
L — i {H, ...} is the Liouville operator and (...) denotes a canonical ensemble 
average. According to the spirit of the CSDM, we try to determine it to the 
lowest order in the form [7,30]: 

A k (w) = 27rA k 5(^ - u k ), (6) 

where the unknown parameters A k and c^ k are to be calculated solving the 



6 



first two moment equations (ME's): 




A k (uj) = -t({SZ^S^}) = 2Nm, (7) 



= 4 £( J ( k ') - J ( k - *))((s£sz v ) + 2 (^,Si k ,)) + 2AM 



(8) 



with m = 



Sj) the magnetization per spin. 



To close this system, we should express all the unknown quantities in terms of 
the SD. One can show[7,30] that the transverse correlation function (5' k h S'I k ), 
appearing on the right-hand side of Eq. (8), can be exactly expressed in terms 
of the SD by the relation: 



On the contrary, the longitudinal correlation function \S£Stk) cannot be 
exactly related to the spectral density A k (u;). So, to close the system of ME's 
(7), (8) we should solve another moment problem by introducing a new SD. The 
simplest way to avoid this difficulty is to resort to the decoupling (S^Si^J 



("Sk )'( 'S'-k) = ^ 2 m 2 5k,o, which means to neglect the correlations between the 
Fourier components of the S?. Of course, the approximation is appropriate 
to describe thermodynamic regimes with finite magnetization as under near 
saturation conditions. In Sec. 4 we will introduce a different procedure which 
allows us to describe also, in a reliable way, regimes with zero or near zero 
magnetization. As final step one must express the magnetization m in terms 





Then, taking into account the one-5 ansatz (6), we have 




(10) 
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of the A k (cj). This is not a simple problem for a classical spin system, also 
when S = 1/2 (for which the exact relation S? = 1/2 — Sj~S~ exists in the 
quantum counterpart). However, as shown in Refs. [7,30] within the spirit of 
the SDM, one can use the expression: 



2 S<2 ~ 2W {S£ Si 

m = -, — '—, (11) 

1 — 252^2 J2k \S^S_ 



which is appropriate for all thermodynamic regimes and reduces, correctly, to 
the near saturation relation m ~ S-(SfS]~) /(2S) = S-£ k (5 k + S: k ) /(2SN 2 ) 
arising from the identity S 2 = (S*f + SfSj with \SfSj\ / S 2 < 1. Then, 
from Eq. (10), as a consequence of Eq. (9) and the ansatz (6), Eq. (11) yields: 



m 



3Tm 
N 



^ 

1 _ Tm V J_ 

1 NS 2 ^ 



;i2) 



With the above ingredients, Eqs. (7)-(8) become a closed system to be solved 
self-consistently. By introducing conveniently the dimensionless variables a = 
m/S, T = T/JS 2 , h = h/JS and u k = Uk/JS, the ME's reduce to: 

gk ^ + ^ (k) + |E " M(k "- ) "" W(k,) . da) 

iV k , u) k > 
a 2 = 1 - 3<T *^ (M) 

AT ^ K Uk 

with 

^ ) (k) = E 1 "?°l k,r - as) 

r I* I 

In the thermodynamic limit N — > 00 one must replace the sum XjcO--)/-^ by 
an integral J 1BZ (...)d d k/ (2n) d in the ci-dimensional k-space and hence Eqs. (13) 
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and (14) can be rewritten as: 



d d k 1 

ji _ ' " J1BZ ri7 N 



U ]i = h + a^ p \k)R(k), (16) 
with 

T r d d k> ^)(k-kQ-^)(kQ 

^ (k) - 1 + a Lz — ^nw(k) — • (18) 

The solution of the problem is complicated and one must consider asymptotic 
regimes for obtaining explicit results or use numerical calculations. Moreover, 
we can note that these equations have physical meaning only if a ^ [30]. 
Hence, they cannot describe the critical behavior of the system or a param- 
agnetic phase in zero external field. As we will see, in Sec. 3 this problem can 
be simply overcome with an appropriate modification of Eq. (18). 

3 Low temperature properties 

We first examine analytically the low-temperature solution of Eqs. (16)-(18). 
The expressions of a and u7k to the first order in the reduced temperature 
T, which allow us to capture some relevant aspects of the low-temperature 
physics of the model (1), are given by: 

oJ k ~ h + ftW(k) - T {lP(h)rt p \k) + l ( f\h, k)} , (19) 

(7~1-T/«(H), (20) 



where 



I " ){h) ~Lz(2irYh + n(p)(ky (21) 

rip) ( T^-f ddk ' ^(k')-^(k-kQ , . 
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At first, we assume h ^ so that no convergency problem for the above inte- 
grals in the k-space occurs. To calculate analytically these integrals, one needs 
an explicit expression of the function fi^(k) which, unfortunately, cannot be 
obtained in terms of elementary functions in the whole 1BZ, for arbitrary 
values of the parameter p > d. Nevertheless, for sufficiently small values of the 
external magnetic field, the dominant contribution to the integrals in Eqs.(19) 
and (20) arises from the low wave-vector excitations as close to the critical 
temperature (see Sec. 4). Then, we can obtain an explicit estimate of a and 
u7 k assuming the dominant behavior of Q^(k) in the 1BZ as k — > 0, provided 
that the coefficients of T in Eqs. (19) and (20) remain finite. Keeping this in 
mind, one can show [16,17,28,29] that for p > d we have for f2^(k) the low- A; 
expansions: 



A d kP~ d + B d k 2 + 0(A; 4 ), p + d + 2 

(23) 

C d k 2 ln{A/k) + 0{k A ), p = d + 2. 



The explicit expressions of the coefficients A d , B d and C d , which depend in a 
cumbersome way on the dimensionality d, the exponent p and a wave-vector 
cut-off A related to the geometrical definition of 1BZ, are inessential at this 
stage and will be omitted. However, for case d = 2, they will be explicitly 
given in Sec. 4 where a comparison with some analytical and Monte Carlo 
predictions is performed. 

Taking into account Eq. (23), the integral I d P \h) can be explicitely estimated 
and we have the following low-temperature representation for the reduced 
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magnetization 



a ~ 1 - T < 



I p f 1 d P ■ 

A 21 V ' P-d' p-d' 



A d AP-'" 



d < p < d + 2 



KdJ^dkk^ 1 h + C d k 2 \n(A/k) ',p = d + 2 



-i 



$ 2^(1,1,1 + p>d + 2 



(24) 



where 2Fi(a,b,c; z) is the hypergeometric function, K d = 2 1 ~ d n~ d / 2 /T (d/2) 
and T (z) is the gamma function. An analogous calculation of lj[\h, k) in the 
expression (19) for cJ k is rather complicated. However, an explicit estimate 
of aJk based on the expansions (23) is irrelevant for next developments. The 
low-temperature susceptibility x can be now easily obtained from Eq. (24) 
by derivation with respect to the reduced magnetic field h. For the reduced 
susceptibility x = x/J we have: 



__ da 



^2^(2,^,^,-^), d<p<d + 2 

K d J A dkk d ~ l \h + C d k 2 \n(A/k)] ~ 2 , p = d + 2 ( 25 ) 

^(2,1,1 + f,-^), p>d + 2. 



The above low-T expressions have a physical meaning for T and h in the near 
saturation regime (a ~ 1) and also for h — > when long-range order occurs. As 
mentioned before, the integrals in Eqs. (19) and (20), and hence the functions 
in Eqs. (24)- (25), could diverge in the limit h — > for particular values of the 
exponent a and the dimensionality d of the lattice. When this is not the case, 
from Eq. (20) a spontaneous magnetization at a finite temperature should 
arise signaling the occurrence of LRO. Of course, the low-k behaviors (23) 
determine the constraints for convergency of the integrals involved in the low- 
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T expressions when h — > 0. It is easy to see that, in the limit h — > 0, the integral 
(21) converges only for <i < p < 2c? with d < 2 and for p > d with d > 2. Then, 
for these values of p and d, a spontaneous magnetization m (T) = o" (T) S 
exists and hence LRO occurs at small but finite temperature, with the result: 

(7 (T)=(7(r,0)~l-T/ d w (0) (26) 

where 1^(0) is a finite quantity whose estimate can be immediately obtained 
from Eq. (24). On the contrary, for a > 2d with d < 2, the integral (21) 
diverges as h — > and no finite solution for a exists at T ^ 0. This means 
that, for these values of a and d, no LRO occurs at finite temperature. The 
global situation is schematically shown in Fig. 1. Here we also distinguished 
the domains where, as we will show in the next Sec. 4, the system exhibits a 
critical behavior like for a Heisenberg model with a nearest-neighbor exchange 
coupling (short-range interaction (SRI) regime) and the long-range nature of 
interactions (LRI) becomes effective. 

In conclusion, our explicit low-temperature results suggest that a transition 
to a FM phase at finite temperature occurs in the domains of the (p — d)- 
plane where LRO exists. In the remaining domains a different scenario takes 
place with absence of a phase transition. These predictions are quite consis- 
tent with the recent extension [17] of the known Mermin- Wagner theorem [19] 
to quantum spin models with LRI's of the type here considered. Other ther- 
modynamic properties can be also obtained within the framework of CSDM 
when the solutions of the ME's are known. Indeed, for our spin model, one 
can easily show that the internal energy u and free energy / per spin are given 
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by [30]: 



„(T,A)=<^ 

= 'J/ {S ° } ~ 2¥P £ J (k) K* 5 -") + (^ S1 -)] • < 27 » 
/(T,fc)=/„ + Ijf 

= A + ^ E ^ jf <w [{sis: k )„ + (s^ k )J . (28) 

Here, if/ is the interaction part of the spin Hamiltonian, /o is the free energy 
per spin for a magnetic model without interactions and (...) j/ denotes a canon- 
ical average as a function of the interaction strength J'. As we see, within our 
approximations all the quantities in Eqs. (27) and (28) can be expressed in 
terms of the transverse SD Ak(cj) and hence all the relevant thermodynamic 
quantities of our classical spin model can be evaluated. In particular, for the 
reduced internal energy u — uj .IS 2 we have: 

t- To r d d k J(k) J(0) 2 

u = -ha / - — -7^^ —a 2 . 29 

J JiBZ (27r) d U k 2J V ; 

Then, in the low temperature limit, from (24) we find: 

ucT-h- 1 -^ + 0(T 2 g(h)) (30) 

where the explicit expression of the term O (T 2 g (jij ) can be obtained in a 
straightforward but tedious way. Hence, for the reduced specific heat Ch = 
(du/dT^j h = C h /S, as expected for a classical spin model [30], we have: 

C h ~l + 0(T 2 g(h)). (31) 
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4 Near zero-magnetization regimes 



As we mentioned before, the decoupling \S£Si k } (S£) ^ s suitable 

for regimes with nonzero magnetization (<r ^ 0). Hence, the basic equations 
used in the previous section do not allow us to explore near-zero magnetiza- 
tion domains in the phase diagram as the critical region. To overcome this 
difficulty, one is forced to find a more appropriate decoupling procedure for 
the longitudinal correlation function which allows us to obtain self-consistent 
ME's appropriate for describing regimes when a — > and preserves also the 
simplicity of the one (^-function ansatz for the transverse SD Ak(a>). A possi- 
ble and successful solution to this problem was suggested several years ago for 
spin models with short-range interactions [7,31]. These studies showed that a 
suitable decoupling procedure when the magnetization approaches to zero (see 
also Ref. [3] for the quantum counterpart) consists in writing (see Eq. (8)): 



As we see, with Eq. (32) only the SD A k (cj) is involved and, inserting in 
the ME (8), one finds for the reduced dispersion relation a; k the same formal 
expression (16) but with R(k) replaced by: 



Notice that, the effect of the decoupling (32) corresponds essentially to achieve 
in Eq. (18) the transformation l/(trJV)X; k , (...) — ► (<r/iV) £k' (•••) [7,30,31]. 
Now, the new ME's can be properly used for an estimate of the main critical 




N 



1 



(32) 




(33) 
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properties of our classical spin model when it exhibits LRO and of the low 
temperature paramagnetic susceptibility in the remaining domains of the (p — 
ci)-plane when no LRO exists. 



4-1 Critical temperature and critical behavior 



In the limit h — > with a > 0, the system of Eqs. (16)-(17) with i?(k) given 
by Eq. (33), become 

cj k = o-tt ip \k)R(k), (34) 
R(k) _ 1 I T [ ^ ) (k-kQ-^)(kQ 

rt w 1 + fi (a)( k ) ( 2vr )rf fi(p)(k')i?(k') ' 1 } 



2 1 - 3TQ(T) 



1 - TQ(T) 



(36) 



where 



Q(T) Lz(2vr) d ^)(k) J R(k)- (37) 
The reduced critical temperature T c (with T c = JS 2 T C ) can be determined 

by imposing the condition a(T c ) = and hence by solving the self-consistent 

equation: 

1 - 3T C Q(T C ) = 0. (38) 

For an explicit estimate of T c , one can again calculate the integrals in Eqs. 
(34)-(37) assuming for f^(k) the dominant contribution as k — > and solving 
our self-consistent equations by iteration. To first level of iteration we find 



l + T c /«(0) 



Then, Eq. (38) yields: 

1 a/rw 



T < = rrwT-; (4°) 
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where 1^(0) is given by Eq. (24). Thus, a critical temperature exists when 
1^(0) is finite and hence in the domains of the (p — <i)-plane where LRO 
takes place (see Fig. 1), as expected. Of course the estimate for the critical 
temperature can be systematically improved by calculating the integral 1^(0) 
using next /c-powers in the expansions (23). 

Here we consider explicitly the two-dimensional case for 2 < p < 4 also for a 
comparison with recent analytical and Monte Carlo results. In this case, the 
coefficients in the expansions (23) are given by [16]: 



A, 



r 2 (p) sin[ir(p-2)/2] 



, 2 < p < 4 



(41) 



5 2 = 2^ 2 c(i-i) [c(i-U)-c(i-f,i 



where ((z, a) = X^o( n + a )~ z ls the generalized Riemann zeta function and 
£(z) = C(z, 0) is the ordinary Riemann zeta function. Then, for the integral 
7^(0), in the case of interest 2 < p < 4, we have: 



4 p) (o) 



A 2 -p T 2 (p) 



sin 



n(p - 2) 



and hence: 



T c (d = 2) 



4-p2 1 -Pir 2 

(4 - p) tt 2 Ap- 2 



2 



(42) 



(43) 



2PT 2 (p) sin 

The critical temperature as a function of p (2 < p < 4) for d = 2 and S — 1/2 
is plotted in Fig. 2 and compared with the corresponding results recently 
obtained for the quantum Heisenberg model by Nakano and Takahashi, us- 
ing a modified spin-wave (SW) theory [14] and the two-time Green function 
EMM within the Tyablikov decoupling [16], and by Vassiliev et al. with a 
Monte Carlo simulation [18]. Our result appears to be consistent with the 
ones obtained for the quantum counterpart in view of the known feature that 
a quantum spin model can be reasonably approximated by a classical one only 
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in the large-S" limit [32]. 



From Eqs. (34)-(37) it is easy to obtain also the behavior of cr(T) as T — > T c + . 
With m = aS and T = JS 2 T, we get (for all values of p and d in the domains 
where a phase transition is allowed): 



m ~ (44) 

J c 



with 



= ^ (45) 



Next, for the reduced paramagnetic susceptibility, defined &sx = (J /h = x/J 
in the limit h — > with a — > 0, Eqs. (16)-(17) with the new expression (33) 
for R(k) reduce to 



_- f d d k 1 

1 = 3xT U(2^rrww (46) 

V ; V ; A 7ibz (2vr) d l + x^(k') V ; 



with X(k) = n^(k)R(k). As usual, an estimate of x(T) as T -> can be 
obtained assuming the low-k behavior (23) for f2^(k) with the aim to find 
a solution of Eqs. (46) and (47) such that x(T) -> oo as T -> T+. Then, we 
find: 



X = Jx~(^t^) (48) 
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with 
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7£r4, U < p < < 



2d~p ' 2 1 



2d, d<2 
d + 2, d > 2 



lin, V = |d , d < 4 



(49) 



1, d <p < < 



\d, d < 4 



d + 2, d > 4 



and 



3*2> 2<d<4 



7=1 lm, d = 4 , P>d + 2. 



(50) 



1, d > 4 

The symbol xi n , here and below, denotes the main (T — T c )-dependence with 
a logarithmic correction (for instance, x ~ — ^c) In 1/ — T c ) ). 

We now determine the behavior of the reduced magnetization a along the crit- 
ical isotherm as h — > 0. Starting from the basic ME's (16)-(17) it is immediate 
to see that, for small values of h and a, the equation for the critical isotherm 
has the form: 

1 - |x 2 - 3T c Q c (a/h) = (51) 

where 

d d k 



Qc (, /5 ) = (, /5)j ( b 



(2^)" 1 + (<r/A)X c (k) (52) 
and X c (k) is determined by Eq. (47) with x replaced by a/h and T = T c . 

Eq. (51) can be solved numerically but an reliable estimate of the reduced 

magnetization a as h — > can be simply obtained, as usual, assuming in 

Eqs. (51) and (52) the dominant contribution of fi^(k) as k — > 0. With this 
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assumption and the definition a ~ h 5 for the isotherm critical exponent 5, 
Eqs. (51) and (52) yield: 



5= I 



2d-p ' 2 



, M < p < < 



2d, d<2 
d + 2, d > 2 



3 ln , p=\d, d<4 



(53) 



3, d <p < < 



|d, d < 4 



rf + 2, d > 4, 



and 



5= <^ 



SI, 2<rf<4 



3ln, 



rf=4 ,P>^+2 



(54) 



3, d > 4 



where 3i n characterizes the behaviour a ~ h 



1/3 



ln/i 



1/3 



Finally, we calculate the critical exponent a for the specific heat in the domains 
of p and d where the transition to a FM phase occurs. It is important to note 
that we cannot use the expression (29) for the reduced internal energy per spin 
which has been obtained assuming the decoupling (S^St^ ~ (S^) (-S'ik)- 
Rather, we must use the general relation (27) with the decoupling procedure 
(32) appropriate near the critical point. Bearing this in mind, the reduced 
internal energy per spin near the critical point assumes the form: 



J(0) 



u ~ -ha - '^^-(T 2 - \to fl + a 2 ) , 

2J 2 V J J 1BZ fa)* cJ k 



d d k J(k)/J 



(55) 



where now Wk involves the quantity R (k) given by Eq. (33). Working for 
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h -> 0, a -> with = x (T) as T -> T*, setting cr ~ x CO h in Eq. (55) 
and taking into account Eq. (46) for x{T) with X (k) evaluated as k^ 0, we 
have for the reduced zero- field specific heat C s=0 (t) = C h=0 (T) / S: 

1 , 1__2 







"2 + 6 X 9T' (T ^ Tc) - 



(56) 



This interesting expression allows us to determine the specific heat critical ex- 
ponent a. Indeed, with x(t)~A(T- T c ) " 7 or \ (t) — A (T - T c ) _1 In [l/(T - T c ) 
Eq. (56) yields: 



C, =0 (T) * 



|-M 2 (t-t c )- q , 



a = 1 - 7 < 



(57) 



i_ IA 2 (T-T C )° In" 1 [5^], a = ln . 



With these definitions and the values of the exponent 7 determined before one 
can immediately obtain the desidered values of a. We find: 



a = < 



2d-p ' 2 U ^ ^ ^ \ 



Oln, P = §<f , 



0, d<p < < 



2d, d<2 
d + 2, d > 2 
d < 4 



§d, d < 4 



d + 2, d > 4, 



(58) 



in the domain where LRI's are active, and 



d A 2<d<4 



a = < 



0l„, 



d=4 , P>d+2 



(59) 



0, d > 4 
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in the SRI regime. One can easily check that the critical exponents (3, a, 7, 
and 5 calculated above satisfy the well known scaling laws 7 = (3(5—1), 
a + 2(3 + 7 = 2 and a + (3(1 + 5) = 2 for all values of p and d which allow 
a transition to the FM phase. Besides, they coincide with those obtained by 
Nakano and Takahashi [16] for the quantum counterpart, consistently with 
the universality hypothesis [32]. 



4-2 Low-temperature paramagnetic susceptibility forp > 2d and d < 2 



The low-temperature behavior of the susceptibility in absence of LRO, i.e. for 
p > 2d and d < 2 (see Fig. 1), is given by Eqs. (46) and (47) but bearing in 
mind that now T c = 0. Following a procedure similar to the one used in the 
preceding section for obtaining the behavior of the paramagnetic susceptibility 
as T — > T c + , in the low-temperature limit, we obtain: 



X ~ S 



exp 



A d K d 
3T 



p — d 
rp~ p-'ld 



T 



T 



p = 2d 

2d<p<2 + d 
(in [^])^~ d p = 2 + d 
p>2 + d. 



(60) 



The previous behaviors generalize those obtained in Ref. [30] for d = 1. It is 
worth noting that for case p = 2d one finds a pure exponential divergence as 
T — > as for one-dimensional model [30]. This result, although consistent with 
that one obtained within the modified spin-wave theory [14] and the EMM 
for the two-time Green functions [16] for the quantum counterpart, does not 
reduce, for d — 1, to the expression derived by Haldane [12] which contains 
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a factor proportional to T . So, one must expect a rf-dependent power law 
factor in T in Eq. (60), which corrects the pure low-temperature exponential 
divergence. The simple lowest-order approximation in the CSDM, here used, is 
not able to capture this important physical aspect. However, within the spirit 
of the SDM [3], one can hope to improve systematically this result working to 
higher order approximations. 

5 Concluding remarks 

In this work we have applied the SDM, within the framework of the classical 
statistical mechanics (CSDM), to study the thermodynamic properties of a 
classical d-dimensional Heisenberg FM model with LRI's decaying as r~ p {p > 
d) in the presence of an external magnetic field. To lowest order in the CSDM, 
the most relevant magnetic quantities have been obtained analytically in the 
low-temperature regime under near saturation conditions as a functions of d 
and p. The FM LRO at finite temperature has been shown to occurs in a wide 
region of the (p — <f)-plane (see Fig. 1) where a transition to a FM phase takes 
place decreasing the temperature. 

The thermodynamic regimes with near zero and zero magnetization have been 
explored on the basis of a proper modification in the ME's of the transverse 
SD. So we have estimated the critical temperature and the main critical prop- 
erties of the model, beyond the mean field approximation, in the domains of 
the (p — d)-plane where a FM ordered phase is expected. The critical expo- 
nents, here determined as a function of p and d, coincide with those ones for the 
corresponding quantum model. This constitutes a microscopic check of the ir- 
relevance of quantum fluctuations, consistent with the universality hypothesis 
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in the theory of critical phenomena. In contrast, as expected, our calculations 
show that the nonuniversal parameters, such as the critical temperature, de- 
pends on the classical or quantum nature of the model under study, especially 
for small values of the spin S as explicitly shown by a comparison of our 
critical temperature estimates with analytical studies and Monte Carlo sim- 
ulations previously achieved for the quantum counterpart. Finally, we have 
determined the low-temperature behavior of the paramagnetic susceptibility 
for p > 2d with d < 2 when no LRO occurs. 

The rich scenario here obtained in a unified and consistent way for a non- 
trivial spin model, shows clearly the potentiality and the effectiveness of the 
SDM also in treating classical many-body systems. It is indeed relevant feature 
that, already to the lowest-order approximation, the CSDM is able to capture 
the essential physics of the model, consistent with Monte Carlo simulations 
[18,28,29] and exact results [20,21,22,23,24], and to obtain results beyond the 
Tyablikov-like approximation for arbitrary values of dimensionality d and the 
decaying exponent p > d. Besides, it offers the possibility to achieve system- 
atically higher-order approximations with the aim to improve the results here 
obtained and to study also the damping of the oscillations on the same footing 
[6]. 

In conclusion, due to the great experience acquired in the quantum many- 
body theory, we believe that the CSDM [5,6,7,8,30] and the formalism of the 
two-time Green functions in classical statistical mechanics [4] constitute a 
promising tool to explore equilibrium and transport properties of a wide va- 
riety classical many-body (not only magnetic) systems. Further developments 
and applications along this directions are desiderable. 
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Fig. 1. Domains of the (p — <i)-plane when a ferromagnetic long-range-order (LRO) 
exists and is absent (NO LRO). The dashed line p = d + 2 separates the domains 
where long-range interaction (LRI)- and short-range interaction (SRI)-regimes oc- 
cur. The dashed region (p < d) corresponds to a nonextensive thermodynamics. 
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Fig. 2. The estimated critical temperature T c /J as a function of the decaying 
exponent p for a two-dimensional spin-1/2 Heisenberg ferromagnet. The full line 
(T c — CSDM) represents the estimate here obtained. The dashed (T c — GF) and 
(T c — SW) lines refer to Green function (GF) method [16] and modified spin-wave 
(SW) theory [14] results. The dots depict the quantum Monte Carlo (T c — QMC) 
predictions [18]. 
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